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Abstract. In thermonuclear supemovae, intermediate mass elements are mostly produced by distri buted burning provide d that a 
deflagration to detonation transition does not set in. Apart from the two-dimensional study by Ropke & Hillebrandt ( 2005), very 
little attention has been payed so far to the correct treatment of this burning regime in numerical simulations. In this article, 
the physics of distributed burning is reviewed from the literature on terrestrial combustion and differences which arise from 
the very small Prandtl numbers encountered in degenerate matter are pointed out. Then it is shown that the level set method 
continues to be applicable beyond the flamelet regime as long as the width of the flame brush does not become smaller than 
the numerical cutoff length. Implementing this constraint with a simple parameterisation of the effect of turbulence onto the 
energy generation rate, the production of intermediate mass elements increases substantially compared to previous simulations, 
in which the burning process was stopped once the mass density dropped below 10 7 gcirr 3 . Although these results depend 
on the chosen numerical resolution, an improvement of the constraints on the the total mass of burning products in the pure 
deflagration scenario can be achieved. 
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1. Introduction 

In the course of the last few years, observational indications in 
favour of a delayed detonation in type la supernovae (SNe la) 
have mounted. For example, calculations of the X-ray spectrum 
of the Tycho supernova remnant assuming various hydrody- 
namical models ap pear to support a d eflagration to detonation 
transition (DDT) (Bad enes et al.l l2006). Furthermore, an inves- 
tigation of nea r-infrared emi s sion li nes of three branch-normal 
supernovae by iMarion et al. I d2006l) implies very little carbon 
residuals at radial velocities less than 1.8 ■ 10 4 cms _1 . Even 
the most advanced three-dimensional simulations of thermonu- 



clear supernovae assuming pure deflagrations ( Ro pke et al 
2006; Schmidt ^: Niemeverl2 006) fail to satisfy this constr aint. 



In models with delayed detonations dGamezo et alj 120051) . on 
the other hand, the supersonic propagation of burning fronts 
dispose of virtually all carbon except for the outermost layers. 
Furthermore, a gravitation al confined detona tion was suggested 
as an alternative scenario (Ple wa et alj|2004l) . 



However, a recent numerical study by Ma ier & Niemevei 



(2006) demonstrated that detonation waves fail to penetrate 
processed material stemming from the initial deflagration 
phase. Therefore, pockets of unburned material ar e likely to 
surviv e even a delayed detonation. Apart from that, Niemeverl 
19991) pointed out several theoretical arguments against DDTs. 



In additio n, accommoda ting the observed variability of SNe 
la dStritzinger et a i1 l2006h within the DDT scenario appears 



to be difficult, because delayed detonations tend to produce at 
least a solar mass of iron group elements and explosion en- 
ergies in excess of 10 51 erg, which are typical values charac- 
terising bright SNe la. But there are also SNe la of moder- 
ate or low luminosity producing much smaller masses of iron- 
group element s and less explosion energy. Apart from that, 
Jha et al. d2006h inferred non-negligible amounts of carbon at 



low expansion velocities from the late-time spectroscopy of 
the SN 2002cx, a peculiar supernova with very low luminos- 
ity. If this result was confirmed by further observations, one 
might conceive of a sub-class of SNe la originating either 
from Chandrasekhar mass white dwarfs which are only par- 
tially burned by pure deflagrations or from sub-Chandrasekhar 
mass progenitors. In this article, we shall be concerned exclu- 
sively with the Chandrasekhar mass scenario. 

Schmidt & Niemever (2006) showed that the explosion en- 



ergy and mass of iron group elements in thermonuclear su- 
pernova simulations with pure deflagration can be varied over 
about an order of magnitude if non-simultaneous point igni- 
tions are applied. To that end, the simple MLT mod e l of the pre- 
supernova core proposed by IWunsch & Woosleyl d2004 was 
adopted for the implementation of a stochastic ignition proce- 
dure. In a numerical case study, models with a total number of 
ignitions ranging from a few up to several hundred events per 
octant were investigated. The main result is that the total mass 
of iron group elements can be adjusted to any value smaller 
than O.75M0 with a maximal explosion energy of roughly 
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0.8 • 10 51 erg. Therefore, these models are feasible candidates 
for less energetic SNe la. 

Due to the artificially chosen termination of thermonu- 
clear burning at mass densities below 10 7 g cirT 3 , however, the 
prediction of the carbon an d oxygen residuals in the m odels 
with stochastic ignition by ISchmidt & Niemevei (2006) was 
not reliable. The transition from the flamelet regime to the 
regime of distributed burning is expected to occur at a mass 
density of about 3 ■ 10 7 g c m 3 (Nie meyer & Woosleylll997 : 
Niemeve r & Kersteinll 19971). Since th e level set method-as im- 
plemented bv lReinecke et aL ( 1999bl) -was applied for the nu- 
merical flame front propagation, the treatment of the distributed 
burning regime remained unclear. As a first approximation, 
distributed burning was mostly suppressed by introducing the 
aforementioned density threshold. This resulted in an overes- 
timate of the amount of unburned material. Apart from that, a 
delayed detonation might be triggered in the late bu rning phase 



1997). However, 



dNiemever & Wooslevlll99 7; Khokhlov et al.l. 
the arguments pu t forward by Niemeverl ( 1999 ) and the nu- 
merical results bv lBell et al.ld2004 shed serious doubt on this 



proposition. For this reason, it appears even more important 
to consider the possibility of the subsonic distributed burning 
mode at lower densities. 

A first attempt to include distribute d burning in deflagratio n 
models of SNe la was made by [Ropke & Hillebrandtl (12005b . 
They carried out two-dimensional simulations demonstrating 
that the explosion energy increases and the fraction of carbon 
and oxygen at low radial velocities is reduced. However, their 
prescription of the turbulent burning speed in the distributed 
regime suffered from a misconception of the relevant scales. In 
this article, we will argue that the existing treatment of burning 
fronts with the level set method can be carried over to the dis- 
tributed regime provided that the burning time scale is smaller 
than the eddy turn-over time scale associated with the numer- 
ical cutoff length. Consequently, this constraint has to be ob- 
served in numerical simulations of thermonuclear supernovae 
with the level set method. 

After reviewing the physics of distributed burning in 
Section|2] we will show in Section[3]that the turbulent burning 
speed continues to be determined by the magnitude of turbulent 
velocity fluctuations at unresolved scales after the onset of dis- 
tributed burning. In Section|H we will discuss results from nu- 
merical simulations of thermonuclear supernovae with stochas- 
tic ignition, where the burning process was terminated based 
on a comparison of the burning and unresolved eddy turn-over 
time scales. Thereby, it was possible to increase the production 
of intermediate mass elements significantly compared to previ- 
ous simulations with a density threshold of 10 7 g cirT 3 for ther- 
monuclear burning. In particular, we found that the total mass 
of intermediate mass elements appears to be independent of the 
details of the ignition process, although slower ignition implies 
a larger amount of fuel at the transition from the flamelet to the 
distributed burning regime. 

2. Distributed burning 

In combustion physics, two different regimes of deflagration 
are distinguished: On the one hand, the flamelet regime, in 



which the microscopic flame propagation speed is solely de- 
termined by the thermal conductivity of the fuel. In the dis- 
tributed burning regime, on the other hand, the transport of 
heat and mass is influenced by turbulence even at scales 
comparable to the w i dth o f the flame brush. As argued by 
iNiemever & Kerstein ( 1997 ). the correct criterion for the tran- 
sition from burning in the flamelet regime to distributed burn- 
ing is that the flame width 8 is about the Gibson length Iq. The 
Gibson length is implicitly defined by the equality of the ve- 
locity v'{{) associated wi th turbulent eddies of size ( and the 
laminar flame speed si am ( Petersf l988): 



v'{( C ) = S\ a 



(1) 



In the flamelet regime, the flame front propagates sufficiently 
fast through eddies smaller than the Gibson length such that 
turbulence has virtually no effect on the internal structure of 
the reaction zone. The width of the flame is then given by 



5 = C s y/xT^ 



(2) 



where x is the thermal conductivity, C& a constant dimension- 
less coefficient and the thermonuclear reaction time scale is de- 
fined by 

ne 

(3) 



Here e is the specific thermonuclear energy release and Q the 
rate of energy generation per unit time and unit volume. The 
time scale r nuc and the conductivity^ also determine the lami- 
nar flame speed: 



X 



(4) 



Once lc ~ 6, however, turbulence will begin to affect the 
transport of heat due to turbulent mixing of preheated mate- 
rial into fuel. As a consequence of the enhanced diffusivity, 
the flame front is broadened. If the turbulence intensity does 
not become too high, however, only the preheat zone will be 
broadened while the reaction zone is not disturbed by tur- 
bulent eddies due to the increased viscosity near the flame. 
Th is was, for instance, dem onstrated in numerical simulations 
by iKim & Menonl llOOOb). The effect of turbulence onto the 
flame can be quantified in terms of a turbulent diffusivity x* 
and the width 6* of the broadened preheating zone. According 
to Damkohler's hypothesis, the propagation speed sj" of the 
broadened flame is enhanced in proportion to the square root 
of the turbulent d iffusivity relative to the thermal diffusivity 
dDamkohlerlll940h : 



J lam 



1/2 



(5) 



This relation follows exactly if the rate of energy generation 
Q is unaffected by turbulence. Since the width of the reaction 
zone tends to become significantly smaller than 6* (by a fac- 
tor greater than 10), this mode of burning is called the thin- 
reaction-zones regime. 

The conductivity x can be related to the viscosity v by 
X = Prv, where the Prandtl number is a dimensionless char- 
acteristic number of the fuel. In stark contrast to most terres- 
trial fluids which have a Prandtl number of the roughly unity, 
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Pr <g 1 for degenerate carb on and oxygen in white dwarfs 
(Nan dkumar & Pethickl 1 1 9841) . Introducing the turbulent vis- 
cosity v* = v t + v, the turbulent diffusivity can be expressed 
as X* - ( v t + V )/P r *> where v t accounts for the effect of turbu- 
lent eddies in the range of length scales from the Kolmogorov 
scale 7/k to the 5*. Thus, equation © can be written in the form 

^ = (|L)'V + 1 )"\ m 

Slam \Pr / V V / 

Assuming that turbulence becomes asymptotically 
isotropic towards scales small compared to the integral 
scale, we may use Pr* ~ 1 for the turbulent Prandtl number 
dYakhot & Orszagl Il986h and the velocity fluctuations v'(i) 
at sufficiently small length sca les I asymp totically obey the 
Kolmogorov-Obukhov 1/3-law ( Frischlfl995 ). 



v'(i) ~ (e{) 1/3 . 



(7) 



Numerical evidence for asymptotic isotropy and the above 
power law in simulations of Rayl eigh-Taylor-u n stable ther- 
monuclear flames was reported by IZingale et alj d2005l) . The 
rate of energy dissipation e = v 3 where the Kolmogorov 
scale 771c is about the size of the smallest eddies. Thus, substi- 
tution of (e(c) 1 ^ on the left hand side of equation ([TJ yields 



.1/3, 



1 

Pr 



4/3 



Since (c ~ 6 upon the onset of distributed burning, it follows 
from the above relation that 6* > 8 ~ Pr" 3/I4 ?7K. In thermonu- 
clear supernovae, 10~ 5 < Pr < 0.1 (Nandk umar & Pethickl 
1984 . Thus, the the broadened flame is much wider than the 
Kolmogorov scale. 

Assuming approximate local equilibrium of the energy flux 
through the cascade of turbulent eddies from size 6* down 
to the Kolmogorov scale, we have v t /v (5*/?7k) 4 ^ 3 » 1. 
Equation (|6]l for the enhanced flame propagation speed then 
implies 

^lam \Pr / \T1k) 

Note that the asymptote si am corresponds to 6* - 

Pr~ 3 / 4 77K which is consistent with the estimate in the previ- 
ous paragraph. Since the definition of the flame width (|2]) im- 
plies x 1 ^ 2 x S, it follows form Damkohler's relation © that 



lam 



Ma 



S* IS. Combining this result with relation (0, it is 
possible to eliminate the width of the broadened flame 6* and 
calculate sJV as a function of known parameters. For a detailed 
specification of the coefficients and the resulting expression for 



lam 



we refer to lKim & Menonl d2000a) 



In terms of time scales, the thin-reaction-zones regime is 
characterised by 



jc_ ji_ 

v'(( G ) ~ Pr 1/2 v' 



(10) 



where the expression on the very right follows from the 1/3- 
law (Q and relation dHJ with 5 ~ £g- For constant viscosity, 
growing energy flux corresponds to a decreasing Kolmogorov 
scale. However, as the time scale Pr^'^^/v becomes much 



smaller than r nuc , turbulent eddies will penetrate and even- 
tually break the reaction zone. In this broken-reaction-zones 
regime, turbulence entrains preheated fuel with patches of al- 
ready burned material at a time scale smaller than the nuclear 
time scale ©. It appears that the dilution of burning material 
within the reaction zone would imply a modification of the rate 
of energy generation per unit volume. As a simple parametri- 
sation, we set Q* = C*Q, where Q is the corresponding energy 
generation rate if turbulence was absent. In the case of pre- 
mixed combustion, which evidently applies to thermonuclear 
fusion, the entrainment of fuel with ash may inhibit but cannot 
boost the burning process, i.e. C* < 1. In any case, the quan- 
titative determination of C* requires three-dimensional direct 
numerical simulations of turbulent burning with a detailed re- 
action network at different densities and varying turbulence in- 
tensity. Nume rical investigati ons pointing in this direction were 
presented bv lBell et alj (120041) . In conclusion, writing the effec- 
tive reaction time scale as r* = r nuc /C*, we have 



> 



p r l/2 y 



(ID 



in the broken-reaction-zones regime. 



(g) 3. Extended level set prescription 



Rein ecke et al.1 (Il999al) introduced the level set method for the 
modelling of burning fronts in thermonuclear supernova sim- 
ulations. The basic idea of this approach is to represent the 
flames by the zero level set of a distance function that is de- 
termined by a partial differential equation. The intrinsic prop- 
agation speed in the flamelet regime is the laminar burning 
speed ii am . Because simulations of SNe la are essentially large 
eddy simulations, however, the effective propagation speed is 
asymptotically proportional to the magnitude of unresolved 
turbulent velocity fluctuations, i.e. t/(A), where A is the res - 
olution of the numerical grid dNiemever & Hillebrandtll 19951) . 
For this reason, the effective propagation speed of the numeri- 
cally computed flame fronts is called the turbulent flame speed 
and denoted by s, . For the calculation of s t , we have a dopted 
Pocheau's model (Poch eadll99HlSchmidt etal]|2006bl) : 



St = «la 



|l+C t |^ 



(12) 



For fully developed turbulence, the flame propagation speed is 
asymptotically given by s t VCt<?sgs- LES of turbulent com- 
bustion in a box indicated that a sound burnin g rate was ob- 
tained with C t not differing too far from unity dSchmidt et al. 
2005). We adopted the value C t = 4/3 which is consistent with 



Peters d!999h . The subgrid scale turbulent velocity q sgs ~ v'(A) 
is comput ed by means of the loc alised subgrid scale model pro- 
posed by dSchmidt et alj |2006a). As a great benefit for the ap- 
plication to thermonuclear supernova simulations, this subgrid 
scale model does not rely on certain flow properties such as 
isotropy or a particular mechanism for the production of turbu- 
lence on large scales. 
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According to lKim & Menonl (I2000al) . the above model for 
the turbulent flame speed can be readily extended into the thin- 
reaction-zones regime: 



' lam 




(13) 



Here, the laminar flame speed is substituted by the enhanced 
flame propagation speed s* introduced in Section (|2]). Since 
s i*am i ncor P ora tes the effects of turbulence at length scales t < 
6* , the contribution of turbulent velocity fluctuations in this 
range of length scales has to be subtracted from q sgs . This re- 
sults in the reduced subgrid scale turbulence velocity <7* gs cor- 
responding to the range of length scales between the width of 
the broadened flame, 6*, and the numerical resolution. In the 
case of Kolmogorov scaling (Q, 



?sgs 



'"'a 



2/3 



1/2 



(14) 



The transition from the flamelet to the distributed burn- 
ing regime is expected to occur within t he range of densities 
from 1 7 to 10 8 g/cm 3 . The calculations of lTimmes & Woosleyl 
d 19921) show that 6 assumes values between 10 2 and a few 
cm within this range of densities. The width 6* of the broad- 
ened flame cannot exceed 6 by more than a fact or ~ 10 witho ut 



turbulent eddies breaking up the reaction zone dKim & Menon 
2000al) . Thus, it appears that 6* < 10 2 cm for the thin- 
reaction-zones regime in thermonuclear supernovae. Even in 
the most elaborate numerical simulations, on the other hand, 
A > 10 5 cm. Thus, q* differs from q sgs by a few per- 
cent at most and the broadened flame may still very well 
be represented by a sharp discontinuity. Moreover, the en- 
hanced flame speed is small compared to q* gs g sgs , be- 
am = 6*16 5, 10 and s lam ~ 1 km/s for 
10' g/cm 3 (Timmes & W oosleyl 1 1992b . whereas q sgs ~ 



cause s'/si 



100 km/s dSchmidt et al.l l2006bl) . Therefore, the flame speed 
model ( fT2l can be carried over unaltered to the thin-reaction- 
zones regime with an accuracy better than about 10 %, which 
is within the intrinsic uncertainties of the model. 

As the exploding white dwarf expands further and the den- 
sity drops below ~ 10 7 g/cm 3 , the nuclear time scale rises 
rapidly and, inevitably, the condition (fTTT) for the broken- 
reaction-zone regimes will be met. The broadened flame will 
then dissolve into a flame brush in which fuel and ash are mixed 
by turbulent eddies. Although there is no well defined width 6* 
of the flame brush, we can specify a turbulent mixing length 
^bum which corresponds to the typical size of eddies with turn- 
over time of the order of the effective reaction time r* : 



kburn 



(15) 



y'(^burn) 

The length scale f\, um can be regarded as generalisation of the 
Gibson length. 

Since the width of the flame in the thin-reaction-zones 
regime is very small compared to typical grid resolutions in 
simulations of thermonuclear supernovae, the flame brush in 
the broken-reaction-zones regime will initially not be resolved 



either, i.e. ^bum ;$ A. In terms of time scales, r* = r nuc /C* < 
A/<7 sgs . As long as this constraint is satisfied, it appears to be 
a sensible approximation to propagate the unresolved turbulent 
flame brush with the level set method. For 



> 



C*A 

Qsgs 



(16) 



on the other hand, the level set representation breaks down. 
Thus, without implementing an explicit treatment of volume- 
burning in the code, the burning process has to be terminated 
in a numerical simulation once the above condition is met. Of 
course, this implies a dependence of the termination of the 
burning process on the numerical resolution. Without includ- 
ing volume-burning, however, this implication is inevitable. 
Former numerical simulations of thermonuclear supernovae 
with the level set method merely appeared to be converged be- 
cause of an artifical termination criterion that was chosen to be 
independent of the resolution (see the following Section). 

The extension of the level set method into the distributed 
burning regime as outlined in this Sec tion markedly differs 
from t he flame speed model applied by iRopke & Hillebrandtl 
(2005), which is based on Damkohler's hypothesis (O and, 
consequently, is restricted to the thin-reaction-zones regime. 
Essentially, their model is based on the assumption s t Sj am 
or, equivalently, 6* - A. It should become clear from the above 
discussion, however, that this assumption does not hold. 

Table 1. Power-law fit of the nuclear energy timescale r nuc in 
the range 0.01 < pg < 0.5. The numerical data are taken from 
Tim mes & Wooslevl d 1992b . 



P9 


Tnuc (numerical) 


T mc (fitted) 


0.5 


0.00089 


0.00126 


0.2 


8.66 • 10~ 6 


3.06 • 10-" 


0.1 


1.18 ■ 10~ 7 


2.29 • 10~ 7 


0.05 


1.76- 10~ 8 


1.71 • 10~ 8 


0.01 


5.23 • irr 10 


5.55 ■ 10- 10 



4. Numerical simulations 

We investigated the effect of including burning with the 
extended level set p rescrip tion in a follow-up study of 
ISchmidt & Niemeverl Q2P06:) . To that end, we modified the 
criterion for the termination of the burning process in the 
code as follows. From the values of s\- dm and 6 calculated by 
iTimmes & Woosleyl d 1992b . we obtained r nuc = 5/siam and lin- 
early fitted log r nuc as a function of logp9 , where p9 is the mass 
density in units of 10 9 g/cm 3 , for py < 1.0. The power law re- 
sulting from this fit is 



c(P9) 



4.146- 10" 

pV av 



(17) 



The data points as well as the corresponding values of the fit 
function are listed in Table Q] Evaluating T nuc (pg) locally in 
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t [s] t [s] 

Fig. 1. Single-octant simulations with N = 128 3 , C e = 10 4 and various choices of C*. Plotted are the evolution of the total energy 
and the integrated mass of intermediate mass elements, respectively. Also shown is a reference simulation, in which the burning 
process was stopped at the density threshold p<) < 0.01. 



Table 2. Total energy release and masses of burning products at t — 5.0 s for the simulations shown in Fig. [TJ 







£ nl , [10 M erg] 


£ km [10 M erg] 


M Nl /M 


M Mg /M 


Pi 


= 0.01 


1.261 


0.753 


0.707 


0.188 


c 


= 0.01 


1.235 


0.727 


0.708 


0.155 


c* 


= 0.1 


1.287 


0.779 


0.708 


0.217 


c* 


= 1.0 


1.332 


0.824 


0.708 


0.270 



each cell, it can be checked whether the criterion [16] is ful- 
filled for a certain prescribed constant C*. If the nuclear time 
scale becomes larger than the threshold given by the right- 
hand side, burning will be terminated. Since A ~ 10 6 cm and 
Isgs ~ 10 7 cms _1 after about one second, C*A/g sgs < 0.1s. 
Thus, the power law [17] has to be somewhat extrapolated to- 
wards pg as 0.003, which appears to be reasonable. 

To begin with, we computed a series of single-octant simu- 
lations on grids consisting of N = 128 3 cells with the code de- 
scribed in lSchmidt et al.l d2006bl) . Due to the hybrid geometry 
of these grids, with a fine-resolved uniform inner part and ex- 
ponentially growing cell size in the outer part, in combin ation 
with the co-expanding grid technique introduced by (Ropke 
20051) . it was possible to study trends even with a relatively 
small number of grid cells. If we had chosen a higher reso- 
lution, the applicability of the level set method would have 
been constrained even further. This can be seen from the cri- 
terion dToT ) with the scaling law g sgs ~ A 1/3 . Since the nuclear 
time scale r nuc gradually increases as the explosion progresses, 
lowering the cutoff length A implies that the termination point 
will be met earlier. On the other hand, a certain resolution 
is mandatory for capturing the large-scale production of tur- 
bulence b y Rayleigh-Taylor instabilities. From the resolution 
study in ISchmidt et al.l ( 2006bh and the influence of doubling 



the resolution of thermonuclea r supernova simulations wit h 
stochastic ignition discussed by ISchmidt & Niemeyei (2006), 
we conclude that 128 3 cells per octant are sufficient for arriving 
at a sensible estimate of the total amount of burning products. 

Applying the stochastic ignition procedure formulated by 



Schmidt & Niemeyer d2006l) . we selected a reference model 
with C e = 10 4 , where the parameter C e adjusts the time scale 
of Poisson processes generating the ignition events in thin ad- 
jacent spherical shells. Depending on the choice of C e , the 
outcome of the explosion, in particular, the release of nuclear 
energy and the total mass of iron group elements is varying 
greatly. For C e = 10 4 , the total number of ignition events per 
octant is roughly 100 and the yield of nu clear energy becomes 
nearl y maximal (see Table 1 and Fig. 1 in lSchmidt & Niemeyei 
( 20061) ). We repeated this simulation with the same initial white 
dwarf model and parameter settings, except for the termination 
of the burning process as explained above, using three different 
values of C*. 

The time evolution of the total amount of intermediate mass 
elements (represented by 24 Mg) plotted in the right panel of 
Fig. Q] clearly shows the influence of the burning termination 
criterion. Assuming C* = 1, the total magnesium mass Mm s at 
t = 5.0 s increases by 44 % compared to the reference simula- 
tion (see Table [2] in which burning ceases for py < = 0.01. 



6 



W. Schmidt: On the applicability of the level set method beyond the flamelet regime 




t [s] t [s] 

Fig. 2. Evolution of the total energy and the integrated mass of intermediate mass elements, respectively, for several full-star 
simulations with N = 256 3 . For each value of the exponentiation parameter C e , which controls the number of ignitions events, 
either p9 < 0.01 or r nuc > A/q sgs (corresponding to C* = 1.0) was set for the termination of the burning process. 



Table 3. Total energy release and masses of burning products at t — 5.0 s for the simulations shown in Fig. [2] 







c e 


En m [10 M erg] 


£ kin [10 M erg] 


M Nl /M 


M Mg /M e 


P ; = 0.01 




10* 


0.985 


0.478 


0.523 


0.201 


P ; = o.oi 


5 


10 3 


1.274 


0.766 


0.715 


0.188 


pi = 0.01 




10 5 


1.165 


0.657 


0.647 


0.184 


C* = 1.0 




10 2 


1.095 


0.587 


0.545 


0.293 


C* = 1.0 


5 


10 3 


1.328 


0.821 


0.701 


0.280 


C* = 1.0 




10 5 


1.246 


0.738 


0.645 


0.286 



One should note that this is only a lower bound, because C* — 1 
indicates that the broken-reaction-zones regime has just been 
entered and volume burning at resolved scales, which cannot 
be treated with the present code, would consume even more 
carbon and oxygen. The mass of iron group elements ( 54 Ni), 
on the other hand, remains unaffected which, of course, re- 
flects nickel being produced entirely within the flamelet regime 
at densities pg > 0.05 while magnesium is stemming largely 
from distributed burning at lower densities. Even for C* = 0.1, 
meaning that the extended level set prescription breaks down 
well within the broken-reaction-zones regime, the yield of in- 
termediate mass elements is still higher than in the case of the 
density threshold p\ = 0.01. If C* = 0.01, on the other hand, 
MMg becomes slightly smaller. Since C* « 1 implies a greatly 
reduced energy generation rate, the physical quenching of the 
flame brush would be nearly reached in this case and, as a con- 
sequence, only little volume burning could possibly follow. In 
this case, we can be fairly sure that the actual amount of burn- 
ing products would be more or less what is seen in the simula- 
tion. 

We also performed full-star simulations with N = 256 3 , for 
which A is initially the same as in the single-octant simulations. 



Varying the exponentiation parameter C e , the burning process 
was terminated once r nllc > A/g sgs , i.e. C* = 1.0. The results in 
comparison to the correspondi ng simulations with the term ina- 
tion criterion py < 0.01 from dSchmidt & Niemeverll2006l) are 
plotted in Fig. [2] Because of the substantially larger mass of 
unburned carbon and oxygen that will be left over at the transi- 
tion from iron group to intermediate mass element production 
if C* is smaller and, consequently, the ignition proceeds slower, 
one would expect Mm s to increase with C e . However, the burn- 
ing statistics listed in Table [3]clearly demonstrates that Mm s is 
almost constant for varying C e . Hence, distributed burning ap- 
pears to consume about the same amount of fuel independent of 
the state of the exploding star at the end of the flamelet regime. 

Fig- ErE show contour plots of the total mass density and 
the partial densities of C+O, Mg and Ni, respectively, in two- 
dimensiona l central sections at time t - 5.0 s. In comparison to 
Fig. 2-4 in ( Schmidt & Niemeyerl2006l) . which show the corre- 
sponding density contours in the simulations with the burning 
process terminated if py < 0.01, one can see that the explosion 
ejecta are substantially enriched in intermediate mass elements. 
This is also illustrated by the plots of the fractional masses as 
function of radial velocity in Fig. |3]-|5] Although the residuals 
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Fig. 3. Total and fractional mass densities in a central section Fig. 4. Total and fractional mass densities in a central section 
and the corresponding probability density functions in radial and the corresponding probability density functions in radial 
velocity space for a full-star simulation with C e = 10 2 at t — velocity space for a full-star simulation with C e = 5 ■ 10 3 at 



5.0 s. 



t = 5.0 s. 



of unburned material at low radial velocities, which have been 
plaguing the deflagration models, are still present in each case, 
it is likely the a significant further reduction would result from 
volume burning beyond the break-down of the level set descrip- 
tion. Nevertheless, a non-negligible amount of carbon and oxy- 
gen at radial velocities larger than about 5000 km/s will be left 
over for any conceivable choice of parameters. 

5. Conclusion 

We have argued that the turbulent flame propagation speed in 
thermonuclear supernova simulations is given by the subgrid 
scale turbulent velocity even in the distributed burning regime. 
Moreover, the level set representation of the flame front re- 
mains valid beyond the flamelet regime provided that the width 
of the flame brush is smaller than grid resolution. In terms of 
time scales, this constraint corresponds the termination crite- 
rion dl6l . 



Very little is known about the interaction between turbu- 
lence and the burning process in the broken-reaction-zones 
regime, in which turbulence is mixing fuel and ash faster than 
the nuclear reactions are progressing. Since the break-down of 
the level set description possibly occurs in this regime, we in- 
troduced the dimensionless parameter C* which specifies the 
reduction of the energy generation rate due to turbulent entrain- 
ment of fuel and ash. The limiting case C* 1 corresponds 
to the thin-reaction-zones regime. In this regime, the flame is 
broadened due to the enhanced diffusivity in the preheating 
zone while the reaction zone is not significantly affected by tur- 
bulent eddies. As turbulence increasingly disturbs the reaction 
zone, C* diminishes. For C « 1, burning will be quenched. 

Setting C* to values in the range from 0.01 to 1, we 
performed several supernova simulations with burning be- 
ing terminated once ( fToT i was fulfilled. In these simulations, 
we applied the stochastic ignition procedure described in 
(Sch midt & Niemeyerl l2006t). For C* > 0.1, a greater frac- 
tion of intermediate mass elements was produced as in refer- 
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Fig. 5. Total and fractional mass densities in a central section 
and the corresponding probability density functions in radial 
velocity space for a full-star simulation with C e = 10 s at t — 
5.0 s. 



ence simulations with termination of the burning process at 
the density threshold pq = 0.01. Remarkably, we found that 
O.3M0 of intermediate mass elements was produced in the case 
C* = 1.0 independent of the rapidity of the ignition process. 
Consequently, the ignition process mainly determines the total 
mass of iron group elements, whereas the production of inter- 
mediate mass elements appears to depend solely on the pro- 
gression of distributed burning. 

Although a substantial increase of the amount of iron group 
elements might result from volume burning at resolved scales, 
which cannot be treated within the present methodology, it 
seems unlikely that the remaining fraction of carbon and oxy- 
gen would be consumed completely. For this reason, if clear 
indications of carbon residuals at least in certain type la su- 
pernovae were not found, delayed detonations would be an un- 
avoidable conclusion. Even in this case, however, studying the 
physics of distributed burning in more detail might very well 
help to clarify the physical mechanism of the putative DDT. 
On the other hand, if burning was completed in the distributed 



mode in some SNe la, small-scale numerical models of burn- 
ing in the broken-reaction-zones regime and explicit treatment 
of volume burning in large-scale supernova simulations would 
be the prerequisites for further progress. 
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